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ABSTRACT 

We present model line profiles of [O n] 23727, [Ne hi] 23869, [O in] 25007, [Fe vn] 26087, [Fe x] 
26374, [O i] 26300, Ha 26563, and [S n] 26731. The profiles presented here illustrate explicitly the pro- 
nounced effects that collisional de-excitation, and that spatial variations in both the ionization parameter 
and cloud column density, have on narrow-line region (NLR) model profiles. The above effects were 
included only qualitatively in a previous analytical treatment by Moore & Cohen. 

By making a direct correspondence between these model profiles and the analytical model profiles of 
Moore & Cohen, and by comparing with the observed profiles presented in a companion paper and also 
with those presented elsewhere in the literature, we strengthen some of the conclusions of Moore & 
Cohen. Most notably, we argue for constant ionization parameter, uniformly accelerated outflow of 
clouds that are individually stratified in ionization, and the interpretation of emission-line width corre- 
lations with ionization potential as a column density effect. For comparison with previous observational 
studies, such as our own in a companion paper, we also calculate profile parameters for some of the 
models, and we present and discuss the resulting line width correlations with critical density ( n cr ) and 
ionization potential (IP). Because the models we favor are those that produce extended profile wings as 
observed in high spectral resolution studies, the line width correlations of our favoured models are of 
particular interest. Line width correlations with n cr and/or IP result only if the width parameter is more 
sensitive to extended profile wings than is the full width at half-maximum (FWHM). Correlations 
between FWHM and n ct and/or IP result only after convolving the model profiles with a broad instru- 
mental profile that simulates the lower spectral resolution used in early observational studies. 

The model in agreement with the greatest number of observational considerations has electron density 
decreasing outward from n e % 10 6 cm -3 to n € ^ 10 1 2 cm -3 and, due to collisional de-excitation effects in 
the lowest velocity clouds, it generates broad flat-topped profile peaks in the lines of lowest critical 
density (e.g., [O n] 23727 and [S n] 226716, 6731). Because the observed profile peaks of both low and 
high critical density lines are often very similar, our favored model requires a contribution to NLR 
emission-line spectra from low-velocity, low-density, and low-ionization gas not included in the model NLR. 
Subject headings: galaxies: nuclei — galaxies: Seyfert — line: profiles 


1. INTRODUCTION 

The models presented in this paper are intended pri- 
marily to complement observations presented in a compan- 
ion paper (Moore, Cohen, & Marcy 1996, hereafter Paper 
I). The models are also suitable for comparison with other 
observational studies, such as those of De Robertis & 
Osterbrock (1984, 1986) and Veilleux (1991a). In this paper, 
we illustrate the pronounced effects that collisional de- 
excitation and spatial variations in both the ionization 
parameter (F) and cloud column density(N col ) throughout 
the narrow-line region (NLR) have on NLR model profiles. 
By comparing the model profiles with those observed, we 
infer the extent to which these effects are likely important in 
the observed profiles presented in Paper I and elsewhere. To 
strengthen further the connection with previous observ- 
ational studies, we also calculate profile parameters for 
some of the models, and we present and discuss the 
resulting line width correlations with critical density ( n cr ) 
and ionization potential (IP) in § 5. 

The first treatment of NLR kinematics that allowed for 
gradients in both the electron density and the cloud velocity 
across the NLR, and which included photoionization calcu- 
lations rather than making very simple assumptions regard- 
ing emission-line luminosity across the NLR, was that of 
Whittle (1982). Since Whittle’s work, there have been two 
important observational results adding new constraints on 


NLR models, suggesting that a new treatment of NLR kine- 
matics is appropriate. 

First, the very high dispersion (10 km s -1 ) and high 
signal-to-noise ratio (S/N) profiles of Veilleux (1991a) show 
well-resolved extended wings. This requires not only line 
emission over a wide range of cloud velocity, but also a 
shallow radial dependence for the line emission so that 
neither the lowest nor the highest radial velocity clouds 
dominate the contribution to the line profile (see Moore & 
Cohen 1994, hereafter MC). Second, we showed in Paper I 
that a wide range in velocity can be inferred from both low 
and high n cr lines and from lines with very different sensi- 
tivity to ionization parameter. Together, these two observ- 
ational results greatly constrain NLR models. In fact, the 
models discussed here are somewhat overconstrained in 
that no single model can satisfy all the observationally 
imposed constraints. Specifically, the model that is consis- 
tent with the greatest number of observational consider- 
ations cannot account for the fact that, in some objects, the 
similarity in profiles of both low and high n CT lines applies to 
the profile cores as well as the wings. This is discussed 
further in §§ 3 and 6. 

2. CALCULATING THE LINE PROFILES 

In this section, we describe our procedure for calculating 
line profiles. In Appendix A of MC, we described this pro- 


301 


302 


MOORE & COHEN 


Vol. 470 


cedure for the purpose of comparison with the analytical 
model profiles presented there. An appropriate portion of 
the description is summarized here for convenient reference. 

A spherically symmetric NLR is assumed in all the fol- 
lowing models, even though biconical emission is clearly 
resolved in some nearby narrow-line regions (NLRs). In 
MC, we showed that the profile wings are not very sensitive 
to either the opening angle or to the orientation of a bicon- 
ical structure containing all the observed line emission, pro- 
vided we rescale the velocity field so that the maximum 
cloud velocity that contributes to the line profile is fixed. 
Thus, as long as we are cautious regarding inferences based 
on the observed profile cores, the results of our spherically 
symmetric models are applicable to a biconical NLR. Little 
information is lost by considering only spherically sym- 
metric models and the corresponding profile wings because 
it is unrealistic to assume purely radial cloud flow for the 
lowest as well as the highest velocity clouds. In this paper, 
we will also demonstrate that the predicted effects of spatial 
variations in the electron density, ionization parameter, and 
column density are sufficiently pronounced that our most 
important conclusions will not likely be invalidated by 
including secondary, nonradial components of the velocity 
field. 

In these models, we incorporate explicit photoionization 
calculations to derive line profiles of [O 11 ] 23727, [Ne hi] 
23869, [O m] >15007, [Fe vn] 26087, [Fe x] 26374, [O i] 
26300, Ha 26563, and [S n] 26731. For each line, the follow- 
ing integral is evaluated numerically : 

hv z )= I 1 <r[n(x), r(x)] 

%) r min/r max 

X ^ 0 [^] obsc x )> x ]* 2 dx ■ W 

In equation (1), s is the cloud size, n is the proton number 
density at the front of the cloud, n cl is the number of clouds 
per unit volume throughout the NLR, v z is the measured 
projected velocity of interest, and t>(x) is the radial velocity. 
The cosine of the polar angle is /i(r z , x), which equals vjv{x), 
as follows from the delta function that selects the projected 
velocity of interest (MC), and 0 = 1 if its argument is less 
than 1 in absolute value and 0 otherwise. Here T is the 
ionization parameter. As in MC, we adopt Ferland’s (1989) 
definition of the dimensionless ionization parameter, T, as 
essentially the ratio of incident ionizing photons to protons. 
Formally, 

r _ jv^y/ftv^V 
4nn H r 2 c 

where v 0 is the threshold for ionizing hydrogen, n H is the 
proton number density (the definition includes both ionized 
and neutral hydrogen), and r is the distance from the pho- 
toionizing source. 

The flux, cr, into each emission line is determined from a 
two-dimensional quadratic spline fit to line strengths com- 
puted using the photoionization code CLOUDY, kindly 
provided by G. Ferland. We assumed stratified clouds for 
all models in the two-dimensional grid (see MC for more 
details). As in Paper I, we use the term “ stratified cloud ” to 
denote a cloud that is individually stratified in ionization. 
We assume s 2 (x) oc n‘ 2/3 (x) (constant mass spheres), n c[ oc 
l/[x 2 t;(x)] (continuity), and that the cloud velocity is 
u(x) x x a . 

The quantity obsc [/x( v ZJ x), x], is a multiplicative factor. 


determined from expressions derived elsewhere (Moore 
1994, 1995), that allows for the effects of both internal and 
external obscuration by dust. To treat internal obscuration, 
we assumed that dust is distributed uniformly throughout a 
cloud and that the line emission is uniform throughout a 
crescent in the interior of the cloud. The size of the crescent 
was calculated using CLOUDY. To treat external obscur- 
ation, we assumed that dust was distributed throughout the 
intercloud medium with a power-law radial dependence. 

In no case were the blue wings of the profiles very sensi- 
tive to details of the obscuration. Therefore, to reduce the 
number of profiles, we restrict the model profiles presented 
here to those in which obscuration is not included. (We do, 
however, include one model in which the treatment of 
obscuration is identical to that of Vrtilek [1985].) The con- 
clusions of this paper are based largely on comparing the 
blue wings of the model profiles to those observed. Because 
obscuration is not included in the models presented here, 
they cannot produce correlations of line shift with IP 
and/or n cr and cannot account for the observed preponder- 
ance of blueward asymmetries. Nevertheless, we have calcu- 
lated profiles for models that include realistic treatments of 
obscuration, and these models do produce blueward asym- 
metric profiles and line shift correlations with IP and n CT . 
The obscuration is mostly internal or external to the emitting 
clouds according to whether the direction of cloud flow is 
inward or outward, respectively. The models we calculated 
in which obscuration was included have some interesting 
implications for the nature of the dust distribution, but they 
are not relevant to the conclusions of this paper (with the 
very important exception, of course, that the blue wings of 
the model profiles are similar to those in which obscuration 
is neglected). In the case of inflow, the blueward asymmetry 
requires stratified clouds. For outflow, the clouds may be 
stratified or or unstratified, even though the obscuration is 
mostly external in either case. (Stratified clouds may be 
optically thick to the incident photoionizing continuum, 
but optically thin to the line emission.) 

We define N col /N min as the column density in units of the 
minimum column density needed for an [O l] 26300- 
emitting zone to exist in the back of the clouds. In some 
models, N col /N min changes in a manner that is inconsistent 
with the assumption, under which <r was calculated, of strat- 
ified clouds throughout the NLR. To make the appropriate 
correction when N col < N min (not completely stratified 
clouds), we calculate the factor by which as 2 (see eq. [1]) 
overestimates a cloud’s contribution to the line profile. We 
formally define r s as the depth at which hydrogen is 50% 
ionized and r fe as the depth at which the [Fe vii] 26087 and 
[Fe x] 26374 lines reach 99% of their final strength. If 
s < r fe , we rescale the Fe line strengths by the factor, s/r fe 
and set all other line strengths to zero. For r fe < s < r s , the 
Fe line strengths are unchanged, [O iii] 25007 and [Ne in] 
23869 are rescaled by the factor, (s — r fe )/r s — r fe ), and the 
[O ii] 23727, [O i] 26300, and [S ii] 26731 line strengths are 
set to zero. Assuming uniform emissivity in the relevant 
zone of the cloud, the emission-line flux from a spherical 
cloud is of the form 

I = n(S- ^ <5 3 ) , (2) 

where S is the depth of the emitting crescent relative to the 
cloud radius (see MC). The preceding scale factors are valid 
to the extent that the second term in the right-hand side of 
equation (2) is negligible. 

Now we present and discuss the grid of detailed kine- 
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matical models. The profile we denote as [O ii] A3727 is 
actually that of [O n] 23726; it is not the blend, 
23726 + 23729. The [O i] 26300 and [S n] 26731 profiles are 
to be interpreted with caution because the flux in these lines 
depends on the cloud column density, which was not cor- 
rected as above for violating the assumption of constant 
mass in the grid of static models. As discussed in MC, the 
fact that emission into the ionized lines changes little as N col 
is increased above N min was exploited by fixing N co] = N min 
+ e, where e < 0.1 dex but is otherwise arbitrary, to reduce 
the dimensionality of the grid of static models from three to 
two (T and n). Although the variation of profile flux with 
decreasing projected velocity in the [O i] 26300 and [S ii] 
26731 model profiles is not accurate in detail, comparison of 
the two profiles illustrates the effects of collisional de- 
excitation because these two lines have very different n cr . [O 
ii] 23727 can illustrate also the effects of collisional de- 
excitation, with the added advantage that the emission from 
a single cloud into this line is not sensitive to the amount of 
gas with N coi > N min . 

Results of models with cloud velocity increasing outward 
are presented in § 3, and results of models with velocity 
increasing inward are presented in § 4. In all the following 
models, we normalize the density profile throughout the 
NLR such that n min = n(r = r max ) = 10 2 cm -3 . This ensures 
that the constraints of low density, implied by [S ii] 226716, 
6731 and [O ii] 23727, can be satisfied. In some of the 
following models, we violate the constraint of high density, 
inferred from the [O in] lines, 24959, 25007, and 24363 (see 
§ 5), by setting n max = 10 4 cm -3 rather than n max > 10 6 
cm -3 . The velocity field is normalized to a maximum value 
of 500 km s" 1 . Neglecting the effects of obscuration, we can 
then parameterize a model NLR by a, y, T(r = r max ), r min , 
and N col {r = r min ), where v oc r a and n oc r y . The parameters 
for each of these models are listed in Table 1. 

Interpreting the profiles is facilitated by comparing with the 
results of the simple analytical treatment of MC. In MC we 
derived the following analytical expression for line profiles: 


I(v z ) = < 


if I v z | > f max , 
if v min < | v z | < m m 
if 0 < | v z | < « min 


where (after normalizing the peak to unity) 


. = i _ M] T i 

vmg l / 

L X^max/ _| L 


(3) 


/? = («' + l)/a, and 

^ core ^wing(^z ^min) ^ * 

The measured projected velocity of interest is v z , and v min 
and r max are the minimum and maximum radial velocity in 
the NLR, respectively. The above assumes a velocity law, 
v oc r®, and an emissivity law, E oc r® , where E is the line 
emission in units of ergs s _1 cm -3 . This expression yields 
profiles with flat-topped profile peaks of half-width u min . 
For the particular case of constant-mass spherical clouds of 
constant ionization parameter obeying the equation of con- 
tinuity, 


The model profiles are presented in Figures 1-13. The 
solid line in each figure represents a model profile after 
convolution with a a = 10 km s ” 1 Gaussian, corresponding 
to the level of thermal broadening of 10 4 K gas, and thus to 
the highest practical spectral resolution of any observ- 
ational study. The dashed lone represents a model profile 
after convolution with a a = 32 km s" 1 Gaussian, corre- 
sponding to the spectral resolution in the red of our own 
study (Paper I). 

3. RESULTS FOR MODELS WITH V INCREASING OUTWARD 

Model L — We begin by presenting a model that repli- 
cates the observed extended profile wings. This model has 
a = 1 (v x r) and y = — 2 (constant ionization parameter). 
The ionization parameter is one throughout this NLR. We 
choose an ionization parameter somewhat higher than the 
value we argue for in MC so that this model corresponds as 
closely as possible with models in which r decreases 
outward (for later comparison [e.g., models 7 and 8]). We 
set r min = 22.5 pc and N c Jr = r min ) = 10 24 cm -2 . The 
column density at r = r min is higher than what is actually 
required to maintain an [O i] zone throughout this NLR in 
order to correspond as closely as possible to a larger model 
NLR in which an outwardly decreasing N col requires a 
higher N col at r = r min (e.g., model 3). There is no obscur- 
ation in this model. Note that n min and F(r = r max ) fix r max at 
225 pc, and n min , y, and r min /r max fix n(r = r min ) = n mai at 10 4 
cm 3 . 

The resulting line profiles are presented in Figure 1. This 
model, which corresponds to an analytical model with 
P = — 1.7 (see MC) and v min = 50 km s _ \ produces profiles 


TABLE 1 

Detailed Kinematical Models at a Glance 


Model 

(pc) 

r 

max 

(pc) 

(cm 3 ) 

OC 

y 

r(r ra J 

T(r„J 

^colfimin) 

(cm -3 ) 

(cm 

1 

22.5 

225 

10 4 

1 

-2 

1 

1 

10 24 

5 x 10 22 

2 

22.5 

225 

10 4 

0.5 

-2 

1 

1 

10 24 

5 x 10 22 

3 

4.1 

411 

to 6 

0.5 

-2 

0.3 

0.3 

2 x 10 24 

4 x 10 21 

4 

4.1 

411 

10 6 

0.5 

-2 

0.3 

0.3 

5 x 10 23 

10 21 

5 

16.4 

411 

6 x 10 4 

0.5 

-2 

0.3 

0.3 

10 23 

10 21 

6 

16.4 

411 

6 x 10 4 

0.5 

-2 

0.3 

0.3 

10 23 

10 21 

7 

71 

710 

10 3 

1 

-1 

1 

0.1 

10 23 

2 x 10 22 

8 

22.5 

2245 

10 4 

1 

-1 

1 

0.01 

10 23 

5 x 10 21 

9 

22.5 

2245 

10 4 

-0.5 

-1 

1 

0.01 

3 x 10 22 

1.5 x 10 21 

10 

2.25 

22450 

10 6 

-0.5 

-1 

1 

10' 4 

3 x 10 22 

7 x 10 19 

11 

4.1 

411 

10 6 

-0.5 

-2 

0.3 

0.3 

3 x 10 24 

7 x 10 21 

12 

4.1 

411 

10 6 

-0.5 

-2 

0.3 

0.3 

10 23 

2 x 10 2 ° 

13 

7.1 

710 

10 5 

-0.5 

-1.5 

1 

0.1 

10 24 

10 22 
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Fig. 1. — The line profiles of model 1, an outflow model with constant ionization parameter. Specific model parameters are provided in Table 1. This 
model replicates the observed extended profile wings but does not include n e > 10 6 cm -3 gas, as required by the [O m] lines in many objects. Solid lines 
represent the model profiles after convolution with a a — 10 km s -1 Gaussian, corresponding to thermal broadening of T — 10 4 K gas. Dashed lines 
represent the model profiles after convolution with a a = 32 km s“ 1 gas, corresponding to the spectral resolution in the red of our own study (Paper I). 


and relative line strengths (other than the [O hi] >14363 to 
[O in] A5007 line strength ratio) that agree with observ- 
ations about as well as any model discussed in this paper. 
However, the range in density and radial coordinates 
encompassed by this model NLR is small. A maximum 
density of 10 4 cm -3 would require an alternative to photo- 


ionization to account for relatively weak [O iii] >15007 emis- 
sion relative to A4363, seen in many objects. Note that [O 
n] A3727 is slightly broader than higher n cr lines because 
collisional de-excitation effects are becoming important in 
the lowest velocity clouds even with n max as low as 10 4 
cm' 3 . 
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Fig. 2.— The line profiles of model 2, an outflow model. This model is identical to model 1 except that the velocity field is less steep. Specific model 
parameters are provided in Table 1. The model includes only a narrow range of radial velocity and profiles not very different from boxcar profiles result. 


Model 2 . — The second model (Fig. 2) is presented to illus- 
trate how a substantial range of radial velocity can be 
inferred directly by visual inspection of the profiles. This 
model is identical to the first except that a = 0.5 (p ac r l/2 ). 
Model 2 corresponds to an analytical model with ft = — 1.3 
and v min = 160 km s _1 . As was evident in the grid of ana- 
lytical models (MC), as emission into the line profiles begins 


to include only a narrow range in cloud velocity, the 
resulting profiles begin to resemble boxcar profiles. 

Model 3— We consider model 3 as (Fig. 3) as the proto- 
type of a model that shows correlations of line width with 
neither ionization potential nor critical density (see § 5). The 
profiles of model 3 are very different from those of model 2 
because model 3 encompasses a much wider range of radial 
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Fig. 3. — The line profiles of model 3, an outflow model. This model includes a wider range of cloud velocities than does model 2, resulting in more 
pronounced wings in the former. Model 3 is otherwise identical to model 2 except for a modest increase in the cloud column density and a modest decrease in 
the ionization parameter, the latter achieved by increasing slightly the distance between the NLR and the central source. Specific model parameters are 
provided in Table 1. We postulate that model 3 is representative of objects that do not show correlations of line width with ionization potential. Note that 
this model cannot account for the similarity in profiles of lines with low and high critical density, observed in some objects. 


velocity. Here v oc r 1/2 and T = 0.3 throughout the NLR. 
The profiles of model 3 are similar to those of model 1 in 
that both models encompass similar ranges of velocity. The 
difference is that the range of velocity in model 1 is sampled 
over a narrow range of radial coordinates and the con- 


straint of high-density gas is violated. 

In model 3, the electron density decreases outward from 
10 6 to 10 2 cm -3 , while N col decreases outward from 
2 x 10 24 cm -2 to 4 x 10 21 cm -2 . Note, however, the effects 
of collisional de-excitation in the core of [O ii] A3727 (as 
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well as of [S n] 26731). This emission line is inconsistent in 
strength and in shape with that observed. If the [O hi] lines 
result from photoionization, which requires the presence of 
n > 10 6 cm 3 gas in at least some objects, then this diffi- 
culty is not easily avoided, however. The constraint of high- 
density gas is discussed in more detail in § 6. 

Note that the wings of [Ne in] >13869 and of [O hi] 25007 
are slightly less extended than those of [Fe vn] 26087 and of 
[Fe x] >16374. This is because the cloud column density 
decreases to N col < N min slightly before the clouds reach 
their maximum velocity at the outer NLR radius. The effect 
can be made either more or less pronounced by varying JV col 
by factors of order unity. 

Model 4— The fourth model (Fig. 4) is the same as the 
third, except that now we have decreased the column 
density of all NLR clouds by a factor of 4, and correlations 
of line width with ionization potential result (§ 5). In this 
model, N c JN m in decreases to less than 1 well before the 
clouds reach their maximum velocity, as is clear from the 
absence of extended profile wings in all but the Fe lines. 
Note that a line width correlation with ionization potential 
would result for a width parameter that is sensitive to the 
wings, but not necessarily for the FWHM. For example, the 
FWHMs of [O hi] >15007 and of [Fe vn] >16087 are nearly 
identical to each other in models 3 and 4, even though the 
full width at zero intensity (FWZI) of >16087 is significantly 
larger than that of >15007 for model 4. 

We postulate that model 4, in which N col decreases 
outward to less than N min well before the clouds reach their 
maximum velocity at the outermost radius, is representative 
of objects that show correlations of line width with ioniza- 
tion potential, and model 3, in which N co , > N min through- 
out the NLR, is representative of objects that do not. Model 
3 can account for the observed profiles of MCG 8-11-11 and 
NGC 4151 discussed in Paper I (no line width correlations), 
while model 4 can account for the profiles of Mrk 79, Mrk 
704, and Mrk 841 (the extent of the wings increases system- 
atically with increasing ionization potential). 

Models 5 and 6 . — Models 5 (Fig. 5) and 6 are presented to 
illustrate how FWHM correlations with Ip/n cr might result 
from low spectral resolution studies but not from high spec- 
tral resolution studies. These models include a somewhat 
lesser range of radial coordinates than do models 3 and 4. 
In addition, we have added a dusty intercloud medium 
(ICM) with an optical depth of 1 across the entire NLR. The 
treatment of obscuration here is identical to that of Vrtilek 
(1985). 

Model 6 (Fig. 6) is identical to model 5, but we have 
convolved the resulting profiles with a Gaussian with 
FWHM = 150 km s 1 (2.5 A at 5000 A) to illustrate the 
effects of decreased spectral resolution. Note that the differ- 
ence in FWHM between [O in] 25007 and [Fe vn] 26087 is 
accentuated by the convolution. The FWHM after convolu- 
tion is 275 km s -1 for 25007 and 331 km s -1 for 26087, 
compared with 250 and 285 km s _1 , respectively, before 
convolution. If we did not know the intrinsic profiles a 
priori and simply corrected for instrumental broadening by 
subtracting the widths in quadrature, we would have 
derived FWHM = 230 km s' 1 for 25007 and 295 km s -1 
for 26087. That is, we would have underestimated the 
FWHM of 25007 and overestimated that of 26087. As in 
model 6, which produces profiles showing a correlation 
between FWHM and IP/n cr (but the wings are no longer 
resolved!), FWHM correlations with IP would thus more 
likely result from low-resolution studies (e.g., De Robertis & 


Osterbrock 1986) than from high-resolution studies (e.g., 
Veilleux 1991a). 

We have noted in § 5.3 of Paper I that Mrk 79 and Mrk 
704 are examples of objects that show a correlation of line 
width with ionization potential if the width is defined by a 
parameter that is more sensitive to extended profile wings 
than is the FWHM. We consider the [O i] and [S n] pro- 
files of Mrk 704 to be "‘narrow” and that of [Fe vn] 26087 
to be “ broad ” (Paper I), but in fact, all three profiles have 
similar FWHMs. We consider the [Fe vn] profile to be 
broad because of its extended blue wing, which significantly 
affects the interpercentile velocity at 20% parameter, IPV 
(20%), used in Paper I. The above applies to models 4 and 5. 
For the NLR corresponding to those models with well- 
resolved profile wings to show correlations of width with 
IP, the width must be defined by a parameter that is sensi- 
tive to extended profile wings rather than by the FWHM. 

Model 7. — Model 7 (Fig. 7) illustrates the difficulty in 
constraining NLR models if we cannot infer a significant 
contribution to the emission lines from high-density gas (e.g., 
from [O hi] and/or [Fe vn] density diagnostics). With- 
out the constraint of both low and high density, we cannot 
rule out an NLR that encompasses only a narrow range 
in radial coordinate, if it nonetheless includes a substantial 
range in radial velocity. For example, such an NLR might 
include only a narrow range of ionization parameters, as 
required by the [Fe vn] 26087 and [Fe x] 26374 profiles, 
even with ionization parameter decreasing outward as 1/r. 

In this model, n oc l/r(r OC 1/r) and v oc r. The cloud elec- 
tron density decreases outward from 10 3 to 10 2 cm -3 , while 
r decreases from 1.0 to 0.1. The ionization parameter 
decreases by only one decade because this NLR encom- 
passes only one decade in radial coordinate. Even this slight 
decrease in T in the direction of increasing velocity 
(outward) has a discernible effect on the [Fe x] 26374 
profile wings. However, the extended wings in all the other 
profiles are roughly consistent with those observed. 

Model 8— Model 8 (Fig. 8) shows the effects of increasing 
the maximum density by even one decade above the pre- 
vious varying ionization parameter model. This corre- 
sponds to an extra decade in the range of radial coordinate 
and, more importantly, in the range of cloud velocity. Here 
t’maxAmin = 100 (it was 10 for model 7), with the result that 
the profile wings are dominated by the low-velocity core. It 
is true in general that if this wide a range of velocity is 
included, the resulting profiles are very unforgiving of any 
deviation from uniformity in the emissivity law, and either 
wing-dominated or core-dominated profiles result (depen- 
ding on whether the emissivity law accentuates the high- 
velocity or low-velocity contribution to the profiles). 

It is important to note that we do not reject this model 
because the profiles are narrow but because almost all the 
flux is contained in the core. If this model NLR were 
observed at realistic S/N, one would infer r max much less 
than 500 km s “ x . If we renormalize the velocity field to, say, 
v max = 2500 km s~ *, so that there is measurable flux at 500 
kms -1 , then r min is increased by a factor of 5 to 25 km s“ b 
The resulting profiles are 5 times broader but represent 
essentially the same model NLR with the property that, 
contrary to observations, the resulting profile wings are 
dominated by the core. 

4. RESULTS FOR MODELS WITH V DECREASING OUTWARD 

Model 9. — Model 9 (Fig. 9) is a conventional gravita- 
tional inflow model (a = —0.5) with density and ionization 
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Fig. 4. — The line profiles of model 4, an outflow model. This model is identical to model 3 except that the cloud column density throughout the NLR is 
decreased by a factor of 4. Specific model parameters are provided in Table 1. We postulate that model 4 is representative of objects that show correlations of 
line width with ionization potential if the width is defined by a parameter that is sensitive to extended wings. 









parameter decreasing outward as 1/r. The [Fe vu] A6087 
profile is reasonable because the ionization parameter at the 
inner radius is 1.0, which is rather high, and even though 
T oc 1/r, T is sufficient to produce significant [Fe vu] >16087 
emission throughout most of this NLR model. The wings of 
this profile are reasonable, in part, because [Fe vu] >16087 
emission from the lowest velocity clouds with the lowest 


ionization parameter is weakened. The core of Ha is strong- 
er relative to the wings than is the core of [Fe vn] A6087. 
Note that the synthetic [Fe x] >16374 line is not at all consis- 
tent with the profiles we have emphasized throughout this 
paper in either shape or strength. In addition, this model 
does not include n 10 6 cm " 3 gas. 

Model 10 . — This model (Fig. 10) illustates the effects of 
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Fig. 5. — The line profiles of model 5, an outflow model. Specific model parameters are provided in Table 1. This model is similar to model 4 but includes a 
smaller range of radial coordinate and velocity, resulting in a correlation of FWHM with ionization potential that is only barely discernible. In addition, due 
to the smaller range in radial coordinate, this model NLR does not include n e > 10 6 cm" 3 gas as required by the [O ill] lines in many objects. We have also 
included obscuration by a dusty intercloud medium. The treatment of obscuration is identical to that of Vrtilek (1985)). 


satisfying the constraints of both low and high density with 
a density profile decreasing as slowly as 1/r The outer 
radius is comparable to that of an entire galaxy and is, 
therefore, completely unrealistic. The flux in the synthetic 
Fe lines is far lower than observations require, and the four 


decades of radial coordinate encompassed by this model 
NLR corresponds once again to v mA Jv min = 100 so that, in 
the other lines, the profile core dominates the wings. 
Although this model produces profiles that are clearly 
inconsistent with those we have emphasized throughout 
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Fig. 6 . — The line profiles of model 6. an outflow model. Model 6 is identical to model 5, but we have convolved the resulting profiles with a Gaussian with 
FWHM = 150 km s~ 1 (2.5 A at 5000 A) to illustrate the effects of lower spectral resolution. The extended wings in the line profiles of high ionization now 
have a discernible effect on the FWHM, as well as on a width parameter that is sensitive to extended profile wings. 


this paper, the profiles may seem not too different from 
those of Mrk 841 (§ 5.4 of Paper I), provided we rescale the 
flux in the Fe lines by a large factor. However, a simple 
scaling of the Fe lines that has little effect on the other lines, 
is inconsistent with the large column density assumed in the 
model. 

Model 11 . — Model 11 (Fig. 11) is a gravitational inflow 


model (a = —0.5) with constant ionization parameter. We 
have already rejected this type of model in the discussion of 
the analytical models (MC) because the flux in the profile 
wings is dominated by that in the core. 

Model 12 . — Model 12 (Fig. 12) is the same as model 11 
except that JV col at the innermost radius has been decreased 
from 3 x 10 24 cm -2 to 10 2 3 cm -2 . Interestingly, this model 
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Fig. 7.— The line profiles of model 7, an outflow model. Specific model parameters are provided in Table 1. Model 7 illustrates the importance of 
constraining independently the range of radial coordinates and/or electron density included in the NLR. We reject this otherwise successful model only on 
the grounds that it violates constraints, inferred from the [O m] lines and the [S n] lines, on the range of electron density encompassed by the NLR. 

yields extended profile wings in [Fe vii] /16087, [Fe x] assumed stratified clouds throughout the NLR. However, 

A6374, and Ha consistent with those observed. This is there is no emission into the other lines from the lower 

because s/r u decreases outward to one well before the velocity-clouds with s < r fe , and profiles resembling 

outermost radius (minimum velocity) is reached, and this boxcars result. 

rapidly decreasing multiplicative factor diminishes the con- Model 13 . — Model 13 (Fig. 13) is a compromise between 

tribution to the profile from clouds of low velocity. This gravitational inflow models that can account for both 

effect was not included in the analytical models, which broad and narrow [O i] ,3.6300 profiles (n xr" 6/5 ; N coi ^ 
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Fig. 8. — The line profiles of model 8, an outflow model. Specific model parameters are provided in Table 1. This model is similar to model 7, but the range 
of density and radial coordinate and velocity is increased. This model illustrates that if the NLR includes a range of radial velocity as large as two decades, 
then it is difficult to produce profiles with comparable flux in the wings and core, as observations require. In model 8, the profile wings are dominated by the 
low- velocity core. 



N min [MC]) and those favorable for strong [Fe vii] 16087 
and [Fe x] 16374 emission throughout the NLR (n oc r~ 2 ). 
In this model, the flux in the wings is weaker relative to the 
core than the profiles we have tried to replicate throughout 
this paper and in MC (with the exception of [O i] 16300 
and [Fe x] 16374), but perhaps not by so much that a 
similar model cannot be reconciled with observations by 


considering small variations in the parameters of this 
model, or perhaps by including secondary effects not 
included thus, far. Similarly, [Fe x] 16734 is weak com- 
pared to [Fe vii] 16087, but perhaps not hopelessly so. In 
addition, although we have argued that NLR gas encom- 
passes a range of four decades in electron density (MC), the 
constraint is not so robust that we can properly rule out a 
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Fig. 9— The line profiles of model 9. Specific model parameters are provided in Table 1. Model 9 is a classical gravitational inflow model with electron 
density and ionization parameter decreasing outward as 1/r, and with velocity decreasing outward as r °' 5 . The profiles of this model are not unreasonable 
except for that of [Fe x] A6374, which is especially sensitive to the decreasing ionization parameter. In addition, this model includes only low-density gas 
(n ral = 10 4 cm' 3 ), contrary to the [O III] line strength ratios in many objects. 


model that encompasses only three decades in electron 
density, as model 13 does. However, the parameters of 
model 13 are very contrived such that even though several 
attributes of the profiles are inconsistent with observations, 
no single attribute is more than marginally so. Only a 
narrow range of parameters can achieve this. 


5. SYNTHETIC LINE WIDTH CORRELATIONS 

In previous observational studies, a common approach 
to inferring the nature of the NLR velocity field has been to 
analyze line width correlations with IP and/or n cr 
(Osterbrock 1981; De Robertis & Osterbrock 1984, 1986; 
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Fig. 10. — The line profiles of model 10, a gravitational inflow model. The parameters for this model are provided in Table 1. Model 10 is similar to model 
9, but model 9 includes a very broad range of radial coordinates so that the constraints of both low- and high-density gas can be satisfied by a density law 
decreasing outward only as 1/r. 


Whittle 1985b; Appenzeller & Ostreicher 1988; De Robertis and the predicted line width correlations resulting from our 
& Shaw 1990; Veilleux 1991b). The success of this approach model profiles have some interesting implications, which we 

depends on the qualitative assumption that lines of high discuss in this section, for the observational studies enumer- 

IP/n cr are emitted more efficiently closer to the central ated above. 

source, relative to lines of low IP/n cr . We have calculated In Figure 14 we plot IPV (20%) versus both n cr and IP for 

explicitly synthetic profiles for a wide range of NLR models, models 3, 4, 9, and 13 (§ 3). Note the considerable scatter. 
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Fig. 1 1 .— ' The line profiles of model 1 1 , a gravitational inflow model with constant ionization parameter. Specific model parameters are provided in Table 
1. We argue against this type of model because the flux in the profile wings is dominated by that in the core. 


even with zero uncertainty in the measured widths. 
Although scatter is expected, even in synthetic correlation 
plots, because the assumption of a specific spread in cloud 
velocities corresponding directly to a specific IP and/or n CT 
was never intended to apply more than qualitatively, the 
amount of scatter is greater than we would have antici- 
pated. We noted in Paper I that our formal errors were well 


within the scatter in the widths calculated from our 
observed profiles. 

As discussed in § 3, we consider model 3 as a prototype of 
a model that shows correlations of line width with neither 
IP nor n CT . Model 4, which differs from model 3 only in that 
the cloud column density is slightly smaller, is the prototype 
of a model that does show correlations of line width with 
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Fig. 12. — The line profiles of model 12. This model is identical to model 1 1 except that the cloud column density is decreased by a factor of 30. Specific 
model parameters are provided in Table 1. 


IP/n cr . However, only if the widths of the lowest n ct lines 
([O ii] A3727 and [S n] >16731) were decreased significantly 
would the line width plot for model 4 (or 3) roughly match 
correlations (or lack thereof) presented in observational 
studies. Here and in Paper I, we have already discussed why 
our favored models require emission from low- velocity and 
low-density gas outside the (parameter space of the) model 
NLR, so the dilemma posed by the large [O n] and [S ii] 


line widths is entirely expected. 

Model 9, a conventional gravitational inflow model with 
n (and T) oc 1/r, yields the tightest correlation of line width 
with IP/n cr of any model we have calculated. Indeed, the 
correlation of IPV (20%) with IP is far tighter than any 
correlation we have seen in the literature. The greatest diffi- 
culty with model 9 is that the total flux in the [Fe] lines is 
far weaker than that observed, desite choosing a high ion- 
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Fig. j 3 _ — The line profiles of model 13, a gravitational inflow model. Specific model parameters are provided in Table 1. This models is a compromise 
model in that several attributes of the synthetic profiles are inconsistent with observations, but no single attribute is more than marginally so. Only a narrow 
range of parameters can achieve this. 


ization parameter specifically to minimize the difficulty. As 
was the case for [O ii] and [S n] in models 3 and 4, the [Fe] 
lines are broad (and weak) because the lowest velocity gas in 
the model does not emit these lines appreciably. Another 
disadvantage of model 9 is that a small change in the 
parameters could not account for an object that does not 
show line width correlations with IP/n cr . Why we invoke 


low-velocity [O ii] and [S n]-emitting gas outside the 
model NLR to reconcile models 3 and 4 with observations 
when we could also invoke low-velocity [Fe vn] and 
[Fe x]-emitting gas outside the model NLR to reconcile 
model 9 with observations is discussed further in § 6. 

In Figure 15 we plot the FWHM of emission lines versus 
both n CT and IP for models 5 and 6. These models are identi- 
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Fig. 14. — The width parameter of Whittle (1985a) (also used and discussed in Paper I), IPV (20%), vs. both n ct and IP for models 3, 4, 9, and 13. The model 
with the strongest line width correlations is model 9, a conventional gravitational inflow model. Our favored models, 3 and 4, yield low n cr line ([O n] >13727 
and [S n] A673 1) profiles that are broader (as well as weaker) than those observed. 


cal except that model 6 has been convolved with a a = 82 
km s” 1 Gaussian to simulate the effects of lower spectral 
resolution. Throughout this paper and Paper I, we have 
considered a profile broad if it shows extended wings. 
However, with high spectral resolution in which the 
extended wings are fully resolved (model 5), the profiles with 
and without extended wings have similar FWHM. Only 
after degrading the spectral resolution (model 6) do corre- 
lations of FWHM with IP/n cr result. 

6. DISCUSSION 

Some of the important results of these models support the 
more qualitative results of the analytical models discussed 


in MC. These results are summarized below. First we 
address some implicit biases in the conclusions inferred 
from these models, and we discuss related observational 
tests. 

6. 1 . Caveats and Related Observational T ests 

An NLR model that encompasses a substantial range of 
velocity over a narrow range of radial coordinate can yield 
profiles similar to those observed (e.g., models 1 and 7). It is 
difficult to constrain such a model. For example, as dis- 
cussed here and in MC, the forbidden Fe lines require a 
narrow range of ionization parameter, but this can be 
achieved even if T oc 1/r, provided r min /r mai is sufficiently 
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Fig. 15. — FWHM vs. both n cr and IP for models 5 and 6. These models are identical except that the profiles of model 6 have been convolved with a a = 82 
km s’ 1 (1.4 A at 5000 A) Gaussian to simulate the effects of lower spectral resolution. The most important effect is that conspicuous FWHM correlations 
result only after degrading the spectral resolution. 


close to unity (see model 7). We can rule out this type of 
model only by invoking the constraints on the range of 
density (MC). Note, for example, that if we increase to four 
decades the range of density included in model 7, then the 
resulting model will include four decades in ionization 
parameters, resulting in Fe line profiles very different from 
those of other lines (see, e.g., model 10). The constraint of 
high-density gas, inferred from the [O in] lines, is not 
applicable to all objects, however. Nevertheless, we assume 
that the same physical processes are at work in objects for 
which we can and cannot infer the presence of high-density 
gas. Therefore, if the constraint of high density leads us to 
conclude that v oc r 1/2 and n e oc \jr 2 in some objects, then 
we assume that these velocity and density laws are applic- 
able in all objects. That is, we assume that objects for which 
we can and cannot infer the presence of high-density gas 
may differ in n max , but do not differ greatly in v(r) and« e (r). 

A related point is that in some objects, the increase in flux 
with decreasing projected velocity is more rapid than for the 
profiles we have emphasized throughout this paper (and 
MC). The most powerful profiles diagnostically are those 
with the most extended wings. When the distinction 


between wings and core is blurred, the observed profiles are 
consistent with a wide range of models. For example, the 
transition from wings to core is more abrupt in the [Fe vii] 
26087 profile of Mrk 704 than of MCG 8-11-11 (Paper I) 
(thus, the distinction between wings and core is defined 
better in the former). The NLR models 5 and 6 are similar 
to model 3 except for a smaller v m ^Jv min . If all NLRs are 
qualitatively similar, then we expect that the same acceler- 
ation mechanisms and basic assumptions apply to all 
objects if they apply to many objects, but differences in the 
line profiles from object to object can be explained, in part, 
by differences in v m Jv mia (or equivalently, r m Jr m J or 
perhaps by small differences in N col (r min ). 

Our favored NLR models (e.g., models 3 and 4) require 
the presence of low-density and low-velocity gas outside the 
model NLR to account for the observed profile cores of 
[O n] 23727 and [S n] 26731. The model profiles of these 
lines differ from those of the other lines because, in the 
superposition of boxcar profiles, the [O ii] and [S n] pro- 
files lack a contribution from the narrowest boxcar profiles 
due to the effects of collisional de-excitation at small r. 
However, the profile similarity observed in lines of both 
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high and low n cr extends into the profile cores in many 
objects (see, e.g., those of Veilleux 1991a). One might rea- 
sonably wonder why we invoke low-velocity [O n] and [S 
ii] emission outside our model to explain the 23727 and 
26731 profiles when we could just as well invoke low- 
velocity [Fe vii] and [Fe x] emission outside the model, 
perhaps relaxing the requirement of roughly constant ion- 
ization parameter throughout the NLR. One reason is that 
if gas outside our model produces low-velocity forbidden Fe 
line emission, this gas must still have an ionization param- 
eter of about 0.3 (assuming photoionization, of course), 
which implies at least some form of continuity with the gas 
producing the high-velocity forbidden Fe line emission. 
[O ii] 23727 and [S n] 26731, on the other hand, are emitted 
over a wide range of ionization parameters. In fact, these 
lines are commonly emitted by gaseous nebulae heated by 
ordinary stars, while [Fe vn] 26087 and [Fe x] 26374 are 
usually associated with active galactic nuclei (AGNs). 

The similarity of the profiles of lines with low and high n cr 
in some objects seems fortuitous if gas outside the NLR 
model accounts for the core of the lowest n cr lines but not 
for the core of the highest n GT lines. However, the difficulty 
posed by the profile similarity is not avoided easily if con- 
straints on the range of electron density encompassed by the 
NLR are taken seriously, and certainly not if the ionization 
parameter is roughly constant. The difficulty is also less 
severe if we note that we have already argued in MC that 
radial cloud flow does not by itself account for the profile 
cores and note also that our favoured model already 
requires some causal connection between the NLR and the 
host galaxy, other than gravitational influence by the latter, 
to explain various correlations of properties of NLR spectra 
with properties of the host galaxy (Busko & Steiner 1988, 
1989; Veilleux 1991c; Whittle 1992a, 1992b, 1992c; Nelson 
& Whittle 1994). We argue for outflow in the NLR, which 
precludes gravity from either the host galaxy or from a 
central black hole as the dominant acceleration mechanism 
of NLR gas. 

Considering that the [O n] 23727 profile core predicted 
by our favored model is inconsistent with observations, 
requiring emission from low-velocity and low-density gas 
outside the model NLR, the validity of the [O in] high- 
density constraint is clearly worth careful attention. One 
important question is whether photoionization of high- 
density gas is the most natural interpretation of relatively 
low values of u = (24959 + 25007)/24363 (such that 25007 
and 24959 emission is suppressed relative to 24363 emission 
by collisional de-excitations). In Moore (1994), it was 
demonstrated that objects with relatively small u tend also 
to have relatively small 25007/H/? line strength ratios, 
arguing (though not strongly) for photoionization of 
high-density gas rather than collisional ionization of low- 
density gas as the more likely interpretation of relatively 
small u. 

A more direct test of the correct interpretation of small u 
is to compare the [Fe vii] line strength ratio, 
(25721 +26087)/25 159, and the [N ii] line strength ratio, 
(6548 + 26584)/25755, with the [O hi] line strength ratio. 
Each ratio is a useful diagnostic of 10 5 <n e < 10 7 cm -3 
photoionized gas and is derived from lines fairly close in 
wavelength. It is unlikely that an alternative to photoion- 
ization could affect the [O in], [N ii], and [Fe vii] density 
diagnostics in the same way. We are proposing a program 
of very high spatial resolution spectrophotometry of spa- 


tially resolvable NLRs on a previously unexplored angular 
scale using the Hubble Space Telescope ( HST ) in the optical 
as well as the Keck 10 m telescope in the mid-infrared. One 
of the goals of the HST observations is to isolate the contri- 
bution to spatially integrated emission-line spectra from the 
highest density clouds, obtaining S/N sufficient to analyze 
the above lines which, though weak, are important diagnos- 
tics of electron density. 

If the cloud electron density is constrained to a range 
such that collisional de-excitational is not important for 
[O ii] 23727, then a constant ionization parameter model 
corresponds to a very narrow range of radial coordinate 
(factor of 10 or so). To sample a range of cloud velocity 
consistent with observed profiles would then require a steep 
velocity gradient; this would rule out gravitational inflow, 
which requires the velocity field to decrease as r -1/2 or 
slower. However, if the NLR has a steep velocity gradient 
and small density range, we do not actually require a con- 
stant ionization parameter. Even with the ionization 
parameter decreasing outward, the [Fe vii] 26087 and 
[Fe x] 26374 profiles are similar to lines much less sensitive 
to ionization parameter because the range of ionization 
parameters is small. If relatively low values of 
u = (24959 + 25007)/24363 can be explained by an alterna- 
tive to photoionization, then we cannot rule out such a 
model (see model 7). It is extremely important, therefore, to 
establish the correct interpretation of relatively small u. 

If gas outside our model produces significant [S ii] 
226716, 6731 emission, then the constraint of low density, 
inferred from the integrated flux in these lines, is not entirely 
valid. In this case, the constraint must be inferred from the 
high-velocity wings in these lines when they exist. Of course, 
extended wings in the [S ii] lines result in blended profiles. 
One of the goals of proposed HST observations is to obtain 
deblended [S ii] 226716, 6731 profiles in an object with 
extended [S ii] wings from spatially resolvable NLRs (e.g., 
NGC 4151, NGC 1068, and Mrk 3). If these observations 
can show that the NLR includes clouds with n e % 10 2 cm -3 
(from the [S n] lines) as well as clouds with n e « 10 6 cm -3 
(from the [O in] and [Fe vii] lines) over a range of two 
decades in radial coordinates, then the observations will 
lend further support to the assumption of constant ioniza- 
tion parameter. In this case, it is reasonable to assume that 
the low-velocity as well as high-velocity forbidden Fe line 
emission is from within the classical NLR even if this 
assumption implies that the low-velocity [S ii] and [O ii] 
emission is from outside the classical NLR. The extra [S n] 
and [O ii] emission might come from H n regions or may 
even come from the outer parts of a dusty torus. The possi- 
bility of low-ionization forbidden-line emission from a torus 
warrants further study. For example, based on observed 21 
cm absorption and on upper limits to molecular absorption, 
it has been calculated that a significant fraction of the gas in 
the dusty obscuring torus in Cygnus A could plausibly be in 
a hot (5000-10,000 K), mostly atomic phase (see Conway & 
Blanco 1995, and references therein). We have assumed also 
that the high-velocity profile wings are emitted within the 
model NLR, whatever the implications for the low-velocity 
profile core. This is a justified bias in that many physical 
processes (some of which are enumerated in MC) certainly 
affect the profiles at the level ofr min . More importantly, the 
observed blue asymmetry constrains the highest velocity 
clouds in a model NLR to move primarily in the radial 
direction regardless of the nature of the velocity field at low 
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velocity. We have also shown (MC) that assuming spherical 
symmetry for a biconical NLR is valid for the profile wings, 
but not for the core. 

6.2. Summary of Conclusions 

We summarize now some of the important results of this 
work. In MC, we argued for constant ionization parameter 
based on the sensitivity of [Fe vn] 26087 and [Fe x] 26374 
emission to the ionization parameter and the similarity of 
these profiles to those of lines much less sensitive to the 
ionization parameter. The modeling of this paper shows 
that if much NLR gas had an ionization parameter outside 
a range of about 0.1 to 1.0, then the Fe line emission would 
differ in its spatial extent from the 25007 emission, for 
example, and, contrary to observations, Fe line profiles very 
different from those of 25007 would result. 

Models 3 and 4 demonstrate that the highest velocity 
clouds must be optically thick to the photoionizing contin- 
uum in objects in which the [O i] 26300 profile has similarly 
extended wings as the [Fe vn] 26087 profile. These two 
models also illustrate how objects that to and do not show 
line width correlations with ionization potential can both 
be explained by a small difference in cloud column density 
of the respective NLRs. In this scenario, N col decreases 
outward as r _4/3 (constant mass spherical clouds with n e oc 
r ~ 2 ) from typically N col % 10 24 cm -2 to N col « 10 22 cm~ 2 . 

Overall, the models presented here support the conclu- 
sion of MC that producing extended profile wings, which 
are often observed, requires that the velocity field increases 
outward if line emission is that of constant-mass spherical 
clouds that are conserved in number and have a constant 
ionization parameter. In our favored models, each radial 
coordinate, r, gives a relative contribution to the line profile 
proportional to r 1/3 which, though not steep, nevertheless 
precludes an arbitrarily wide range of cloud velocity from 
contributing to the model profiles if they are to be consis- 
tent with those observed. For a realistic range in radial 


coordinates of two decades or more, the velocity field must 
then increase outward as about r° 5 rather than about r 1 0 
so that comparable flux is emitted into the wings and the 
cores of the model profiles. 

One can easily imagine physical processes that could 
invalidate our assumptions of cloud continuity and of con- 
stant cloud mass, especially since we argue for outflow so 
that the clouds and ICM need not be comoving. Regardless 
of these assumptions, however, what is required to produce 
profiles similar to those observed can be described simplisti- 
cally as follows. The observed profiles cannot be replicated 
by the boxcar profile corresponding to a single velocity. A 
superposition of boxcar profiles, corresponding to a range 
of velocity, is required. The superposition must be such that 
neither the wings nor the core dominate the profile. That is, 
the flux in the wings is comparable to that in the core. A 
strictly uniform emissivity law can satisfy this condition 
even if the range of velocity is infinite (v min = 0). However, if 
a substantial range of velocity is included, the resulting pro- 
files are very unforgiving of any deviation from uniformity 
in the emissivity law, and either wing-dominated or core- 
dominated model profiles result (depending on whether the 
emissivity law accentuates the high-velocity or the low- 
velocity contribution to the profiles). If the assumption of 
cloud continuity is altered drastically, then another assump- 
tion, such as that of constant cloud mass, must be modified 
in a contrived fashion to compensate such that the resulting 
emissivity law still favors neither the highest velocity clouds 
nor the lowest velocity clouds by a large factor. 
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Note added in proof . — Since submission of these two papers, we have observed NGC 4151 with long-slit spectroscopy on 
the Keck 10 m telescope with 0'.'65 FWHM spatial resolution, comparable to the spatial resolution needed to test some of the 
predictions of both of these papers. For that object, the following predictions of these papers are confirmed: 

1. The [Fe vii] 26087 and [O t] 26300 emission-line fluxes have comparable angular extents on the subarcsecond angular 
scale that partially resolves the source of most of the line emission. 

2. Emission into the high-velocity wings of the [S it] 26716, 6731 profiles coincides spatially with that into the high-velocity 
wings of the [Fe vn] and [O i] profiles. 

3. On the subarcsecond scale, the profiles of all three lines show a velocity shift that increases systematically with angular 
displacement from the nucleus, relative to the systemic velocity. 




